Coulomb Interaction-induced Checkerboard Patterns in Disordered Cuprates 
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We study the effect of the Coulomb interaction on the local density of states (LDOS) and its 
Fourier component in disordered cuprates. It is shown that the Coulomb interaction suppresses 
strongly the maximum value of the LDOS induced by the dopant impurity at each energy and 
expands significantly the Friedel oscillation in real space. The existence of the Coulomb interac- 
tion with a moderate strength yields an energy-dependent checkerboard LDOS modulation around 
the impurity, which is very different from that produced by pure quasiparticle interference. The 
orientation and transformation of the checkerboard pattern with energy and the relations among 
the modulation vectors, dopings and the bias voltages agree qualitatively with the recent STM 
experiments. 
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PACS numbers: 74.62.Dh, 74.72.-h, 74.25.-q 

It is well-known that at half filling cuprates are the 
Mott insulators, where the on-site Coulomb interaction is 
larger than the kinetic energy of electron. When oxygens 
(holes) are doped into cuprates, the electrons can move in 
the Cu02 plane and cuprates become strange metals. In 
other words, the "effective" on-site Coulomb interaction 
decreases. With further increasing the concentration of 
holes up to about 5%, cuprates become superconducting 
and have a d-wave symmetry gap. So cuprate super- 
conductors inevitably possess disorder due to the oxygen 
doping. It is out of question that the electronic states in 
the Cu02 plane are mainly determined by disorder and 
the Coulomb interaction. 

Recently a series of the STM experiments have re- 
vealed the electronic states in both real and momentum 
spaces for the superconductor F^S^CaC^Os+a [1,2,3]. 
It was found that in real space the LDOS has a weak 
energy-dependent checkerboard pattern. At low bias 
voltages, the charge modulation with a period |^ ori- 
ents along 45° to the Cu-0 bonds. At high bias voltages 
well below the superconducting gap, the charge modu- 
lation with a period |j is parallel to the Cu-0 bonds. 
We note that qA are the wave vectors connecting the 
end points of the same Banana-like equal-energy con- 
tours while qs are the wave vectors connecting the near- 
est end points of the neighbor Banana-like equal-energy 
contours [1,2,4,5]. These two kinds of charge modula- 
tions coexist at intermediate voltages. With increasing 
the bias voltage or doping, qA becomes shorter while qg 
becomes longer [1,2]. In Ref. [3], McElroy et al. studied 
in detail the the nanoscale electronic disorder induced by 
the oxygen dopant atoms. It was indeed verified that 
these dopant defects generate the LDOS modulations 
and also suppress superconducting coherence peaks in 
Bi 2 Sr2CaCu20 8+( 5. 

A lot of theoretical investigations have devoted to these 
fascinatingly experimental observations [4-13]. In mo- 
mentum space, the relation among the modulation wave 
vectors, dopings and the bias voltages seems to be well 



explained by quasiparticle interference due to weak disor- 
der. However, in real space, the LDOS near the impurity 
has a stripe-like structure along both 45° to the Cu-0 
bonds and the Cu-0 bonds rather than the checkerboard 
pattern at each energy [9]. In addition, the Friedel os- 
cillation only exists in a short distance (~10 lattice con- 
stant) away from the impurity site. For many random 
impurities, a checkerboard pattern and its transforma- 
tion with energy and doping are observed theoretically. 
But the checkerboard pattern only covers a small inde- 
fective area and is different from that seen by the STM 
experiments, which expands into the bulk. Therefore, it 
is difficult to interpret simultaneously the experimental 
phenomena in both real and momentum spaces by us- 
ing pure quasiparticle interference induced by disorder. 
What factor makes the LDOS at different energy have 
the two-dimensional period structure in the bulk? It is 
natural to consider the Coulomb interaction in the Cu02 
plane as a candidate because the undoped cuprates are 
the Mott insulators. In this paper, we investigate the 
influence of the Coulomb interaction on the Friedel oscil- 
lation in disordered cuprate superconductors. We shall 
see that the experimental phenomena are the results of 
disorder correlation mediated by the Coulomb interac- 
tion. 

The Hamiltonian describing the scattering of quasi- 
particles from a single defect with local modifications of 
both hopping and pairing parameters and the Coulomb 
interaction in a d-wave superconductor can be written as 
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where < i,j > denotes that the sites i and j are 
the nearest neighbor sites, fi is the chemical poten- 
tial to be determined by doping, e k — ti(cosk x + 
cosk y )/2 + t2COsk x cosk y + t^(cos2k x + cos2k y )/2 + 
t4(cos2k x cos k y + cosk x cos2k y )/2 + t$ cos 2k x cos 2k y 
with ii_ 5 =-0.5951, 0.1636, -0.0519, -0.1117, 0.0510 (eV). 
The band parameters are taken from those of Norman et 
al. [14] for Bi 2 Sr 2 CaCu208+5, and the lattice constant 
a is set as o=l. The order parameter away from the im- 
purity is given by A k = Aofcosfe,, — cosk y )/2. Stij = St 
for i = or j = and for others. SAij = SA\ for i = 
or j = 0, SA 2 for i = 1 or j = 1, and for others. U s 
and V m are the nonmagnetic part and the magnetic part 
of an on-site potental, respectively. U is the strength of 
the on-site Coulomb interaction while V is that between 
the nearest neighbor sites. 

The Hamiltonian (1) cannot be diagonalized exactly 
due to the presence of the Coulomb interaction. Because 
the charge modulations observed by the STM experi- 
ments are weak, here we are only interested in the case 
of weak to moderate U and small V, where the charge 
density wave (CDW) and spin density wave (SDW) can 
be neglected [15]. In fact, the presence of CDW and 
SDW only produces an energy-independent peak at the 
modulation wave vectors in the Fourier component of the 
LDOS and doesnot affect the results in this paper [16-22]. 

Following the standard Green's function method [23], 
we can solve the Hamiltonian (1) in the absence of CDW 
and SDW. Then the LDOS measured by the STM exper- 
iments reads 
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where v — 0, 1, N is the number of sites in 
the lattice, and CQ^,{iuj n ) is the Green's func- 
tion in terms of the Matsubara frequencies and 
has the symmetry G kk ,(±ia;„) = G k , k (±iu; n ). 
To the first order in V s , St, SAi and SA 2 , the 
Green's function G? kk ,(iu;„) = G ki ,(iw„)c>kk'<W + 
G kl/ (iw n )G kV ,(iw„){5 kk ^ k [(-l)^'^/k, + £w+i<?k,] 
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(A k + 7k)^k], A k = a k (l - d k ) + bkCk, 
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ously, the weak magnetic potential V m doesnot appear 
in G kk / (iuj n ) because it has only the higher order 
contributions to the Green's function. 

Here we point out that the Coulomb interaction itself 
doesnot produce the charge modulations in the present 
case. However, it has a strong impact on the LDOS in- 
duced by disorder. Substituting the Green's function 
G kk ,(±iw„) into the formula (2), we can obtain the 
LDOS images at different energy and doping. In our cal- 
culations, we have chosen N = 400 x 400, [i = -0.1246 
corresponding to about 15% doping, Ao = 44 meV and 
V 8 = 2St = -2SAx = -48A 2 for simplicity. Here we'd 
like to mention that the other parameter options, includ- 
ing SAi and 5A 2 to be positive and different doping, 
arrive at the similar conclusions. 

Figure 1 shows the LDOS patterns with a moderate 
U = 0.25 eV and a small V = 0.01 cV at different en- 
ergy due to a single impurity at the center of a 40 x 40 
square. In this case, the ratio of the on-site Coulomb 
interaction to the nearest neighbor hopping y « 1.7. 
Obviously, the LDOS patterns for u> > have similarly 
periodic structures with those for u < 0. At u = ±12 
meV, the LDOS has a checkerboard modulation with a 
period ~ 5a. Its orientation is along 45° to the Cu-0 
bonds. With increasing energy to ±18 meV, the modula- 
tion period becomes shorter to ~ 4a. We note that there 
also exists an about 20 x 20 checkerboard modulation 
along the Cu-0 bonds around the impurity at ±18 meV. 
At high energy, e.g. u = ±30 meV, the LDOS shows a 
checkeroard pattern along the Cu-0 bonds, which pos- 
sesses a multi-period structure. Obviously, there doesnot 
exist the long-wavelength modulations. When uo = ±24 
meV, two kinds of the checkerboard patterns along both 
45° to the Cu-0 bonds and the Cu-0 bonds coexist. The 
charge modulation and its transformation with energy 
agree qualitatively the STM experiments. Obviously, the 
Coulomb interaction suppresses the maximum LDOS on 
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FIG. 1: (Color online) The LDOS p(r,u>) at different energy 
due to a single impurity at the center of a 40 x 40 square. 
U = 0.25 eV and V = 0.01 eV. 



or near the impurity sites and magnifies the Friedel os- 
cillation in a large area at each energy [9]. 

In order to further understand the influence of the 
Coulomb interaction on the electronic states in the Cu02 
plane, we study the Fourier component of the LDOS, 
i.e. pq(uj) = Sr e ~ iq 'V( r ! a; )- I n Fig. 2, we present 
Pq(u>) along (jr, 0) direction at different energy. Obvi- 
ously, Pq(ui) has a peak at the modulation wave vector 
qA, which changes with u>, U and V. With increasing U 
and V (u> fixed) or uj (U and V fixed), the peak moves 
to the origin and \qA\ becomes shorter. If U or V has 
an enough large value, the peak will disappear at high 
energy below the superconducting gap. This seems to be 
consistent with the STM experiments [1,2]. In contrast, 
the peak at qa produced by the pure quasiparticle inter- 
ference never disappears when energy approaches to the 
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FIG. 2: (Color online) The Fourier component p q (w) of the 
LDOS along the antinodal direction at different energy and 
Coulomb interaction. 



superconducting gap. 

Figure 3 shows Pq(u>) along (n, tt) direction at different 
energy. We can see that when U <~ 0.1 eV and V = 
eV, the peak at the modulation wave vector qs does not 
move with U at each energy. However, qb becomes longer 
if U >~ 0.1 eV or a small V is presented. Therefore, 
when U and V are fixed, the peaks at qa and qs have an 
opposite shift with increasing energy. In other words, qA 
becomes shorter while qs becomes longer. The relation 
between ^(to) and u> has been observed by the STM 
experiments [1,2]. Note that both qA and qs are very 
sensitive to the presence of V while a small U term has 
a large impact on qA- 

In summary, we have investigated in detail the effect 
of the Coulomb interaction on sites and between the 
nearest neighbor sites on the electronic states in disor- 
dered cuprate superconductors. It is shown that the 
Coulomb interaction plays an important role in form- 
ing the checkerboard patterns of the LDOS modulations 
induced by the impurity. Our model (1) produces the es- 
sential features of the energy-dependent LDOS modula- 
tions as observed in the STM experiments [1,2,3]. In fact, 
cuprate superconductors contain a lot of random oxygen- 
dopant impurities. However, all the impurities correlate 
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FIG. 3: (Color online) The Fourier component /Oq(w) of 
the LDOS along the nodal direction at different energy and 
Coulomb interaction. 

each other through the Coulomb interaction. We con- 
clude that it is nothing but this kind of correlation that 
leads to the energy-dependent checkerboard patterns of 
the LDOS modulations in cuprate superconductors. 
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